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The construction of initial data for black-hole binaries usually involves the choice of free pa- 
rameters that define the spins of the black holes and essentially the eccentricity of the orbit. 
Such parameters must be chosen carefully to yield initial data with the desired physical proper- 
ties. In this paper, we examine these choices in detail for the quasiequilibrium method coupled to 
apparent-horizon/quasiequilibrium boundary conditions. First, we compare two independent cri- 
teria for choosing the orbital frequency, the "Komar-mass condition" and the "effective-potential 
method," and find excellent agreement. Second, we implement quasi-local measures of the spin of 
the individual holes, calibrate these with corotating binaries, and revisit the construction of non- 
spinning black hole binaries. Higher-order effects, beyond those considered in earlier work, turn out 
to be important. Without those, supposedly non-spinning black holes have appreciable quasi- local 
spin; furthermore, the Komar-mass condition and effective potential method agree only when these 
higher-order effects are taken into account. We compute a new sequence of quasi-circular orbits for 
non-spinning black-hole binaries, and determine the innermost stable circular orbit of this sequence. 

PACS numbers: 04.20.-q, 04.25. Dm, 04.70.Bw, 97.80.-d 



I. INTRODUCTION 

Recently, significant progress has been made in nu- 
merically evolving black-hole binaries [1-4]. A major 
goal of these simulations is to estimate the gravitational 
waveform produced by astrophysical black-hole binaries. 
These waveforms will ultimately be used to aid in the 
detection and interpretation of the gravitational wave 
signals we expect to see in observatories such as LIGO, 
VIRGO, TAMA, and GEO600. In order for such simula- 
tions to yield astrophysically relevant results, the initial 
data must be constructed to be astrophysically realistic. 

A very effective approach for constructing numerical 
black-hole binary initial data has been developed and ex- 
plored by two of the authors [5] (see also Refs. [6-8]). In 
Ref. [5], the authors focused attention on two specific, 
limiting cases of binary initial data: corotating black 
holes and irrotational (non-spinning) black holes. These 
numerical initial-data solutions were compared against 
previous numerical results[8, 9] and to analytic post- 
Newtonian estimates [10, 11] for binaries in circular or- 
bits and appear to give the best agreement yet between 
numerical and analytic models of close black-hole bina- 
ries in circular orbits. However, various aspects of the 
physical content of these initial-data sets have not been 
fully tested. 
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In order to construct quasi-circular orbits (as opposed 
to general elliptical orbits), Gourgoulhon et al.[7] pro- 
posed that the data must satisfy a simple condition: the 
Komar[12] and the ADM[13] masses must agree if the or- 
bits are quasi-circular. The Komar mass is only a reason- 
able definition of the total mass of a system if the system 
is stationary. The ADM mass, is an invariant measure 
of the total mass of a system as measured at space-like 
infinity. So, for binary systems that are quasi-stationary 
when they are in quasi-circular orbit, the ansatz seems 
quite reasonable. This ansatz has been tested by com- 
paring numerical models with analytic post-Newtonian 
models. It has also been tested in the case of neutron- 
star binarics[14] where an independent method of deter- 
mining circular orbits exists. For black-hole binaries, an 
independent method of determining circular orbits exists 
in the so-called effective potential method^]. In this pa- 
per, we will further test the Komar- ADM mass ansatz by 
constructing circular orbits using the effective-potential 
method. 

Another important aspect of the binary initial data 
constructed in Ref. [5] that has not been adequately ver- 
ified relates to the spins of the individual black holes. In 
Ref. [5] , the spin of each black hole is fixed by a particu- 
lar choice of boundary conditions applied at the surface 
of the black hole. For the case of corotating black holes, 
the choice of boundary conditions is unambiguous. How- 
ever, for non-spinning black holes, the boundary condi- 
tions were chosen in a way that should be correct in the 
limit of large separation between the black holes in the 
binary. However, this non-spinning ansatz has not yet 
been checked. 

Of course, we must keep in mind that the angular mo- 
mentum (spin) of an individual black hole in a close bi- 
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nary system is not rigorously defined in general relativ- 
ity. To measure the spin, we will have to rely on a quasi- 
local definition. There are many such definitions [15]. For 
our purposes, we will use a definition first made rigor- 
ous by Brown and York[16] and also derived within the 
more recent isolated horizons framework of Ashtekar and 
Krishnan[17, 18]. We will explore extensively the spins 
of the individual black holes in our binary initial data us- 
ing this quasi-local definition of the angular momentum. 
A major result of this study is that we must refine our 
method for setting the boundary conditions in order to 
construct models where the black holes are not rotating. 

We begin in Sec. II with a review of the entire for- 
malism used to construct our initial data. In Sec. Ill 
we discuss various issues associated with computing the 
quasi-local spin of a black hole. In Sec. IV we examine 
the case of corotating binaries, focusing first on exploring 
the criteria for defining circular orbits and then examin- 
ing the spins of corotating black holes. In Sec. V we turn 
to the case of non-spinning binaries and define the correct 
approach for obtaining non-spinning black holes. Finally, 
in Sec. VI we present results related to the inner-most sta- 
ble circular orbits for both non-spinning and corotating 
binaries, and review the major points from the paper. 



II. INITIAL-DATA FORMALISM 

The black-hole initial-data sets that we consider be- 
low are constructed using the conformal thin-sandwich 
decomposition]!!), 20], a set of boundary conditions im- 
posed on the black-hole excision surfaces[5, 6] and at 
asymptotic infinity, and a set of assumptions for various 
freely-specifiable fields. Below, we will outline the most 
important details of the various pieces of our approach. 



A. The conformal thin-sandwich decomposition 

In this work, we will use the standard 3+1 decompo- 
sition with the interval written as 

ds 2 = -a 2 dt 2 + j ij (dx i + /3 i dt)(dx j + ftdt), (1) 

where 7^ is the 3-metric induced on a t — const, spatial 
hypersurface, a is the lapse function, and (5 l is the shift 
vector. The extrinsic curvature of the spatial slice, K^, 
is defined by 



and 
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where C n denotes the Lie derivative along the unit nor- 
mal to the spatial slice, n^. Einstein's equations, in vac- 
uum, then reduce to four sets of equations. Two are 
evolution equations for the spatial metric and extrinsic 
curvature: 



dtjij = -2aKij + 2V ( i/3,-), 



(3) 
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+ P l V e K lJ +2K e{l \7 j) p i . (4) 
The remaining two are the constraint equations 
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Here, Vi, Rij, and R are, respectively, the covariant 
derivative, Ricci tensor, and Ricci scalar associated with 
the spatial metric 7^ . Finally, the trace of the extrinsic 
curvature is denoted K = K\. 

The conformal thin-sandwich decomposition employs a 
York-Lichnerowicz conformal decomposition of the met- 
ric and various other quantities [2 1-23]. The conformal 
factor, ip, is defined via 



7ij = ^lij; 



(7) 



where 7^ is a "conformal metric". The time derivative 
of the conformal metric is introduced by the definitions 



J ij Uij = 0. 



(8) 
(9) 



From this, it follows that the tracefree extrinsic curvature 
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Av = K ij - W 3 K takes the form 
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where a = ip~ e a is the conformal lapse function, and 
u y = UkiJ lk J : ' 1 ■ Furthermore, (LV) is the conformal- 
Killing (or longitudinal) operator acting on a vector, de- 
fined by 



(LI/),, 



(ifj) - §7yVkV , 



(11) 



where is the covariant derivative compatible with 7^ . 
Finally, the conformal tracefree extrinsic curvature can 
be written as 



2<5 V Hl 



(12) 



In terms of our conformally decomposed variables, the 
Hamiltonian constraint (5) can be written 

VV - - J2^ 5 K 2 + l^AijAv = 0, (13) 

where R is the Ricci scalar associated with 7^, and the 
momentum constraint (6) as 



irV i K-%(±u i i)=0. (14) 
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Furthermore, the trace of Eq. (4) can be written as 
V 2 (</> 7 <5) - 0A 7 5) Ur+JL^K 2 + |^- 8 I y l« 



= [d t K - /? fc V fc A-J . (15) 

Within the conformal thin-sandwich formalism, the 
fundamental variables are: ip, a, (3 1 , 7y, K, dt^ij = 5y, 
and c^i-sT. Of these, 7^, A, jiy, and <9 t if represent the 
eight gauge and dynamical degrees of freedom of the grav- 
itational field. These fields must be chosen based on the 
physics of the initial data one wishes to model. The re- 
maining fields, if), a, and f3 l represent the constrained de- 
grees of freedom. Once the other fields have been fixed, 
these fields are determined by solving Eqs. (13-15) as a 
set of coupled elliptic equations. 

Formulating a well-posed elliptic system requires that 
we impose boundary conditions. Typically, these systems 
are solved under the assumption that the spacetime is 
asymptotically flat. If we let r denote a coordinate radius 
measured from the location of the center of energy of the 
system, then as r — ► 00 we have that 



- a\ 



= 1, 

= (n x r y 
= 1. 



(16a) 
(16b) 
(16c) 



f2o is the orbital angular velocity of a binary system, 
or the rotational angular velocity of a single object, as 
measured at infinity. The boundary condition on the shift 
is chosen so that the time coordinate, f = cm' 1 + (3^, is 
helical and tracks the rotation of the system[5-7, 24, 25]. 
If we wish to consider systems with one or more black 
holes, and if we excise the interior of the black hole to 
avoid difficulties with singularities, then we must also 
impose boundary conditions on the excision surfaces. 



B. Black-hole excision boundary conditions 

In this paper, we are interested in the situation in 
which each black hole is in quasiequilibrium and the 
boundary conditions required to achieve this were worked 
out in Refs. [6] and [5]. The assumptions are essentially 
the same as those required of an "isolated horizon" (cf. 
[17, 26-28]). To ensure that the black hole is in quasiequi- 
librium, we enforce the following conditions. First, we 
demand that the expansion 9, of the outgoing null rays, 
k^, vanish on the excision surface, 5, thus forcing the 
boundary to be an apparent horizon: 



= 0. 



(17) 



Next, we require that the shear of the outgoing null 
rays also vanish on the excision boundary, 



a \iv\ s = 0- 



(18) 



In the absence of matter on S, Eqs. (17) and (18) are 
sufficient to imply that 



C k 0\ 



0. 



(19) 



That is, initially, the apparent horizon will evolve along 

Conditions (17), (18) and (19) are coordinate indepen- 
dent, however the final demand breaks precisely this co- 
ordinate freedom. This final condition is that the coor- 
dinate location of the apparent horizon docs not move 
initially in an evolution of the initial data. Let the ex- 
cision boundary surface, S, be a spacclike 2-surface with 
topology S 2 and define s l to be the outward pointing unit 
vector normal to the surface. If we let 



W EE -^(„/* + a M) 



(20) 



represent the set of outgoing null rays associated with S, 
and with the time vector written as 



then this condition can be expressed as 



= 0. 



This immediately yields 



?8i 



(21) 



(22) 



(23) 



To further analyze Eq. (23), we split the shift vector into 
its component normal to the surface, P ± , and a vector 
tangent to the surface, (3/, defined by 



(24) 
(25) 



where hij = Jij — SiSj is the metric induced on S by 
7y. We see that Eq. (23) is a condition on the normal 
component of the shift, 



(26) 



The conformal transformation on 7^, Eq. (7), induces 
a natural conformal weighting for hij and for the unit 
normal to 5, 



1 -2 ~i 
1p S . 



(27) 
(28) 



In terms of these conformal quantities, the condition in 
Eq. (17) takes the form of a nonlinear Robin- type bound- 
ary condition on the conformal factor ip, 
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We now turn our attention to Eq. (18), which can be 
rewritten as 



Dp 



(30) 

where D l denotes the covariant derivative compatible 
with hij. Below, we will assume Ujj = 0; in that case, 
Eq. (30) implies that f3^ 1 must be a conformal Killing 
vector of the conformal metric, hij on the boundary S: 



£ (i /V } 



0. 



(31) 



In practice, we write the parallel components of the shift 
as 



(32) 



where £ l is a rotational conformal Killing vector on 
with affine length of 27r, and Sl r is a constant. As shown 
in Ref. [5], the choice of /3 M * directly parameterizes the 
spin of the associated black hole. The restriction on (3^ 
is quite remarkable. Regardless of the choice of the con- 
formal metric 7^ (and thus for any hij ), Eq. (31) still 
allows sufficient freedom to allow for the parameteriza- 
tion of a rotation about any direction (by the choice of 
£ l ) and with any magnitude (by the choice of O r )[5]. 

To summarize, the quasicquilibrium conditions defined 
in Eqs. (17), (18), and (22) define boundary conditions 
on the conformal factor, ip, via Eq. (29) and on the shift 
vector, /3 l , via Eqs. (26) and (32). These total to four 
of the five necessary boundary conditions for solving the 
coupled elliptic equations associated with the conformal 
thin-sandwich equations. Missing is a condition on the 
conformal lapse, a. 

However, as clearly shown in Ref. [5], the excision 
boundary condition on the lapse is intimately associated 
with a degeneracy in the choice of the initial slicing con- 
dition. In fact, for stationary black-hole spacetimes, it 
is the choice of the lapse on the excision boundary that 
uniquely fixes a particular initial slice. For initial data 
representing systems that are nearly stationary, that is 
systems in quasiequilibrium, it has also been shown that 
the choice of the excision boundary condition for the 
lapse is largely irrelevant. Following Ref. [5], we will 
choose the excision boundary condition for the lapse from 
a set of convenient and rather generic conditions. Assum- 
ing that the excision boundary is spherical, we use 



dai}) 



dr 
daip 



dr 
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(33a) 
(33b) 
(33c) 



C. Quasi-circular orbits 

The formalism reviewed so far in Sees. II A and II B 
provides a very general framework for constructing initial 
data for black-hole binaries. Let us now make a specific 
choice for part of the freely-specifiable initial data. Con- 
sider the choice = d t "fij = 0. This choice implies that 
the conformal three-geometry is, at least momentarily, 
stationary with respect to the time vector, t, given by 
Eq. (21). Our choice for the boundary condition on the 
shift at infinity given in Eq. (16b) implies that our time 
vector has a helical form where the amount of "twist" in 
the vector is parameterized by = |^2o|- 

If we assume that t is an approximate Killing vector 
of the spacetime, then = is a direct consequence. 
With a proper choice of Q , the resulting initial data 
represents a binary system where the black holes are in 
quasi-circular orbits. This statement begs two questions: 
1) How do we determine the proper choice of f2o so as to 
obtain quasi-circular orbits? 2) How do we interpret the 
initial data for other choices of f^o? 

An answer to the first question was proposed by Gour- 
goulhon et al.[7] where they made the ansatz that flo 
be chosen so that the ADM mass[13] and the Komar 
mass[12] agree. We call this the "Komar-mass ansatz", 
and when it is applied, we refer to it as the "Komar-mass 
condition" . 

The ADM energy is an invariant definition of the total 
energy of a spacetime as measured by an inertial observer 
at spacelike infinity. The ADM energy is written as a 
surface integral at infinity 



1 

loTT 



V,- [Qi-m )d 2 s 



(34) 



and fij 



is a flat metric. 



where (/', , = jij - fij, Q = ~i %1 Gv 

If the total linear momentum of the system (as measured 
by the same inertial observer at infinity) vanishes, then 
the ADM energy is usually referred to as the ADM mass. 
In our notation, the Komar mass can be written as 



m k = -L 

47T 



{V l a-^K l ,)d 2 S\ 



(35) 



The Komar mass is a valid expression for the total mass 
of a system only when the system possesses a global time- 
like Killing vector so that the system is stationary. It 
is therefore quite reasonable to assume that by choos- 
ing f2o so that the Komar-mass ansatz will yield initial 
data which is nearly stationary (i.e. in quasiequilibrium), 
which would require that the black holes be in quasi- 
circular orbits. 

The effectiveness of the Komar-mass ansatz has been 
tested numerically in black-hole initial data by com- 
parison with post-Newtonian data for binaries in cir- 
cular orbits[5, 10] (and for neutron-star binaries[14]). 
There are also additional theoretical reasons for expect- 
ing that configurations satisfying the Komar-mass con- 
dition (-Eadm = M K ) will represent systems in quasicqui- 
librium (cf. Refs. [7, 29] and references therein). These 
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arguments show that a system in quasiequilibrium nec- 
essarily satisfies the Komar mass ansatz. However, as far 
as we know, there is no guarantee that a system satisfy- 
ing the Komar-mass condition is necessarily in quasiequi- 
librium. It would therefore be interesting to compare 
the Komar-mass ansatz against an independent method 
for determining circular orbits. The effective-potential 
method [9] will be used below to provide such a compari- 
son. 

It is worth noting that when using the Komar-mass 
ansatz in the context of binary systems, care must be 
used in evaluating Eq. (35). An observer moving along 
the approximate helical Killing vector t is not an iner- 
tial observer. For a true stationary spacetime, the helical 
Killing vector can be split globally into separate timelike 
and rotational Killing vectors and the timelike Killing 
vector is used to define the timeline of the inertial ob- 
server making the measurements. The approximate heli- 
cal Killing vector cannot be split globally, however it can 
be split asymptotically at spacelikc infinity. In evaluating 
the Eq. (35), it is necessary to remove the l~2o x r term 
from the shift, so that the same inertial observer is 
used to evaluate both M K and E ADM . Often, the Komar 
mass is written in a from similar to that of Eq. (35), but 
with the term involving the shift absent. In many cases, 
this yield a correct expression for the Komar mass since 
the contraction of the shift with the extrinsic curvature 
falls off faster than 1/r 2 when the flo x r contribution 
to the shift is omitted. However, it is not always cor- 
rect to simply drop this term. For example, in Painleve- 
Gullstrand coordinates [30-32] and their extension to the 
full Kerr-Newman spacetime [33], this term contains the 
entire contribution to the Komar mass. 

In order to answer the second question of how to in- 
terpret the initial data when Qq is no longer chosen via 
the Komar-mass condition, we must no longer think of 
the helical time vector as an approximate Killing vec- 
tor of the spacetime. A more general interpretation 
of the choice Uij = dtjij = is that the "velocity" 
of the conformal three-geometry vanishes on the initial- 
data slice. In the context of a binary configuration this 
suggests that the system is at either periccnter or apo- 
center of some general bound or unbound orbit. When 
f7 = 0, the system will have no orbital angular momen- 
tum and will represent a generalized version of Misner[34] 
or Lindquist[35, 36] initial data. 



D. Corotating and non-spinning black- hole binaries 

In constructing black-hole binary initial data, we will 
certainly need the ability to specify the spins of the in- 
dividual black holes. In Sec. II B, we mentioned that the 
excision boundary condition on the parallel components 
of the shift, f3 n l , can be used to set the spin of each black 
hole. It is tempting to want to interpret f2 r in Eq. (32) 
as the rotational angular velocity of the black hole. How- 
ever this is not the case. To understand the role of (3/ 



in determining the spin on the black holes, it is useful to 
consider the special cases of corotating and non-spinning 
black holes. While neither case is expected to be seen 
astrophysically, these cases are useful because they repre- 
sent situations where we know either what the boundary 
condition should be or what the final spin should be. 

Corotating black hole binaries represent the case where 
the black holes are rotating synchronously with the or- 
bital motion. A great deal of attention has been payed 
to such binaries because they represent the only config- 
uration of two black holes that can possess a true helical 
Killing field[29, 37-41]. A serious fault with such space- 
times is that they cannot be asymptotically flat since 
they contain a balancing amount of incoming and out- 
going radiation for all time. However, for our purposes, 
they are ideal. First, because we at most have only an 
approximate helical Killing vector, we can still have an 
asymptotically flat solution. More importantly, we know 
that the corotating case possess a Killing horizon. This 
means that the proper choice for the parallel components 
of the shift is unambiguously given by [3^ = 0. 

Since corotating black holes necessarily have some non- 
vanishing rotation (and therefore spin) as seen by an 
asymptotic inertial observer, clearly f2 r cannot represent 
the rotational angular velocity of the black hole. At 
least for a Newtonian binary system in corotation, the 
rotational angular velocity should equal the orbital an- 
gular velocity of the binary system. For post-Newtonian 
computations of corotating black-hole binaries, this is, 
in fact, the condition used to set the spins of the black 
holes[10, 42]. To leading order, this is correct. However, 
as we will demonstrate below, higher order corrections 
are needed to correctly estimate the rotational angular 
velocity of corotating black holes. 

Non-spinning black-hole binaries are, in some sense, 
the simplest case. Here, the individual black holes have 
no spin as measured by an asymptotic inertial observer. 
In previous works, we have referred to such systems as 
having "irrotational" black holes [5, 6]. Perhaps, however, 
this is not the best terminology as it carries with it the 
notion of fluid motion which is not appropriate for black 
holes. From now on, we will refer to such black holes 
as being "non-spinning" . How do we choose /3y % to yield 
black holes with no spin? We argued previously [5] that 
for a binary with orbital angular velocity f^o, we should 
choose tt r — rJo- Also, for spherical excision boundaries, 
we chose £j as the flat space rotational Killing vector, 
projected into the excision surface, that generates a ro- 
tation about an axis parallel to the orbital angular mo- 
mentum vector. We also showed that this condition lead 
to reasonable results. To leading order this is correct, but 
again we will show below that it is necessary to modify 
this choice for the boundary condition in order to pro- 
duce non-rotating black-hole binaries. 
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E. Conformally flat maximally sliced models 

So far, we have discussed the approach used to con- 
struct binary initial data in rather general terms. Wc 
have discussed the conformal thin-sandwich decompo- 
sition and the choice of boundary conditions for con- 
strained variables within that formalism. These con- 
strained variables (ip, (3 l , a) can only be determined after 
values for the unconstrained variables (7^, 5y, K, dtK) 
have been chosen. In Sec. II C, we discussed the inter- 
pretation of initial data when the choice d t "fij = Uij = 
was made. We will use this choice for «y in all models 
constructed below. 

We will also make the assumption of conformal flat- 
ness in all of the models we construct. That is, we will 
choose jij to be a flat metric. This choice impacts upon 
the physical content of the initial data that we will con- 
struct. It is well known that the spatial metric for a 
relativistic binary system cannot be conformally flat [43], 
nor can the spatial metric of a stationary spinning black 
holc[44, 45] or even the metric of a boosted black hole. 
However, the errors introduced by the assumption of con- 
formal flatness are not "grave" . There are efforts under- 
way to improve the choice of the conformal metric. Some 
of these efforts involve using an analytic metric obtain 
from post-Newtonian theory (cf. Rcf. [46]). However 
a more self-consistent approach has been developed by 
Shibata ct al.[47, 48] for the case of neutron-star bina- 
ries and it should be possible to adapt this approach to 
black-hole binaries in the future. 

Finally, we choose to use maximal slicing, K = dtK = 
0, in all of the models we construct. Which member 
of the family of possible maximal slices we choose will 
depend upon our choice for the boundary condition for 
the lapse [5]. As a slicing condition, we do not expect this 
choice to impact significantly upon the physical content 
of the initial data. 



F. Numerical code 

We solve the elliptic equations of the conformal thin- 
sandwich decomposition using the pseudo-spectral collo- 
cation method described in Refs. [5, 49]. As usual with 
spectral methods (see, e.g. [50]), the solution is expressed 
as a truncated series of basis- functions, and is represented 
by a set of expansion coefficients. For appropriate basis- 
functions, discretization errors decay exponentially with 
the number of retained basis-functions [50]. The elliptic 
equations take the form of a set of nonlinear algebraic 
equations for the expansion coefficients. These algebraic 
equations are solved with a Newton-Raphson method, 
where in each step a linear problem is solved via stan- 
dard Krylov subspacc techniques [51] like preconditioned 
fGMRES [52]. 

For the binary black hole solutions of this paper, we 
need to solve elliptic equations in a computational do- 
main with two excised spheres. To do so, the computa- 



tional domain is split into smaller subdomains, namely 
spherical shells and rectangular blocks. A shell is placed 
around each excision surface and a third shell, using a 
compactified radial coordinate, extends from some inter- 
mediate radius surrounding both holes out to the outer 
boundary which is typically placed at a radius of ~ 10 8 . 
The space between is filled by a collection of rectangular 
blocks. 

Further details of the numerical code can be found in 
Ref. [5]. In particular, Fig. 6 of Ref. [5] displays the 
convergence of the code for a typical configuration. The 
calculations below were performed at a resolution com- 
parable to N = 60 in this figure and correspond to dis- 
cretization errors on the order of 10 -5 or 10 -6 for most 
of the quantities we consider. 



III. COMPUTING THE QUASI-LOCAL SPIN 
OF A BLACK HOLE 

In the previous discussion, we referred to the spin of an 
individual black hole in a binary system. However, there 
is no unique, rigorous definition of the spin of a black 
hole unless it is in isolation and stationary. We can. how- 
ever, rigorously define the total angular momentum of an 
initial-data slice as measured by an incrtial observer at 
infinity. Usually referred to as the ADM angular momen- 
tum, it can be expressed as 

J (£) = -L j {Ka - 7ij K) & d 2 S\ (36) 

where £ l is a Killing vector of 7^ . For asymptotically 
flat data, we can choose this to be each of the three flat- 
space rotational Killing vectors in order to determine the 
three components of the total angular momentum of the 
system. 

Although there is no unique, rigorous definition for the 
spin of an individual black hole in general, we can esti- 
mate the spin based upon a quasi-local definition. There 
are many different quasi- local definitions for spin [15]. 
While motivated in different ways, they tend to take 
on a similar form within a 3 + 1 framework. One of 
the earliest useful definitions was derived by Brown and 
York[16]. This same definition was rederived more re- 
cently within the isolated and dynamical horizons frame- 
work of Ashtekar and Krishnan[17, 18]. In cither case, 
the quasi-local spin can be expressed as 

%) = ^ £ i K n - 1*J K ) ? ( 37 ) 

where now is a Killing vector of hy. We note the 
remarkable similarity between Eqs. (36) and (37). The 
main difference between them is that the former is eval- 
uated at infinity while the latter is evaluated on the ap- 
parent horizon. For the initial data we are constructing, 
this will be the excision surface. In terms of the variables 
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used in the conformal thin-sandwich decomposition, we 
evaluate the spin as 



%) = 



1 

Ifhr 



1 



(38) 



where we have used the fact that £ l is tangent to the 
excision surface, i.e. Jij£ l s J = 0. 

The only choice that must be made in evaluating the 
spin is to choose the Killing vector The problem is 
that, in general, an exact Killing vector will not exist. 
However, there are two reasonable choices that can be 
made. If our excision surface (apparent horizon) is a co- 
ordinate sphere within our flat conformal geometry, then 
we can choose to approximate £ J by one of the 3 flat- 
space rotational Killing vectors centered on the excision 
sphere and projected onto its surface. We will denote 
these three choices by 





./•V' r <. 


(39a) 




= x j s e iyj , 


(39b) 


Clz) 


= xW zj , 


(39c) 



where we assume Cartesian coordinates and x s is mea- 
sured relative to the center of the excision sphere. 

An alternative is to attempt to solve Killing's equa- 
tion. Once we have solved the constraint equations and 
have full initial data, we know h%j, the physical metric 
projected onto the excision surface S. In Ref. [27], Drcycr 
ct al. outlined a general method for finding the Killing 
vectors on a closed 2-surface. We give the details of our 
implementation of this method in Appendix A. Here, we 
simply note that exact solutions of Killing's equation will 
not exist in general so what we find are "approximate" 
Killing vectors. 

It is difficult to make a meaningful quantitative mea- 
sure of how far a given solution deviates from being a 
true Killing vector. For our purposes, we will attempt 
to gauge the accuracy of our measured spins by com- 
parison with expected values. For the case of corotating 
black holes, we know that our boundary condition on (3^ 
is correct and we know, at least to leading order, what 
the spin of a corotating black hole should be in a binary. 
Thus, we can compare our quasi-local definition for the 
spin of a black hole against an analytic result. 

Below, when we discuss the computed quasi-local spin, 
we will use the following notation for simplicity: 



Sy 

Sr. 

S K 



computed using Eq. (39a) 
computed using Eq. (39b) 
computed using Eq. (39c) 

computed using an approximate 
solution of Killing's equation. 



(40) 



important test case because it is the one configuration for 
black-hole binaries that is compatible with a true helical 
Killing vector[29, 37, 38]. The thermodynamic relations 
obtained by Friedman et al. [29] , 



SE A 



Sln<5J A 



KiSAi 



(41) 



should apply to our conformally-flat data if we had a 
true helical Killing vector and we assume no local change 
in entropy. In Eq. (41), J ADM is the magnitude of the 
total ADM angular momentum of the system, Kj and 
Ai are respectively the surface gravity and area of the 
Killing horizon of each black hole. Corotating black holes 
also have a non-vanishing spin and, more importantly, 
a physically well-defined notion of the rate of rotation. 
Together, these give us a firm analytic foundation against 
which we can test our initial data. 

Let us first consider the thermodynamic identity in 
Eq. (41). We are free to define a fundamental length 
scale for each of our initial-data solutions. We can use 
this freedom to scale our solutions in an attempt to have 
it satisfy this identity. Of course, Eq. (41) allows for too 
much variation for a single length rescaling to guarantee 
the enforcement of the identity in general. However, we 
can use the freedom to define a fundamental length scale 
x(s) along a sequence of initial data to enforce 



5E, 



fl 5J A 



(42) 



but 



This approach has been discussed previously[5, 
for completeness, we cover it again here. 

Using our freedom to define the fundamental length 
scale x( s ) along a sequence we define the dimensionful 
total energy E ADM (s), total angular momentum J AD m(s)j 
and orbital angular velocity f2o(s) consistently via, 



E AD m(s) = x(s)e(s), 
Jadu(s) = X 2 (s)j(s), 



(43) 
(44) 
(45) 



Then, to enforce Eq. (42), it is sufficient to determine 
the change in x( s ) between two points on the sequence. 
This can be done by integrating along the sequence from 
a point si to another point S2- Doing so, we find 1 



XM = X(si)exp^ - 



! * e>(s)-u(s)j'(s) 
e(s) ~ 2u(s)j(s) 



ds 



(46) 



where a prime denotes differentiation along the sequence. 

Now, if the Komar-mass ansatz is reasonable, then a 
sequence of initial data with varying separation that sat- 
isfies this ansatz should represent a binary in a nearly 
adiabatic evolutionary sequence of quasi-circular orbits 



IV. COROTATING BINARIES 



Although corotation is not considered to be an astro- 
physically realistic state for black-hole binaries, it is an 



1 Note that the similar equation in Ref. [5] contains a factor of 2 
error in the denominator of the integrand. 
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that satisfy Eq. (41). But, having scaled the data to sat- 
isfy Eq. (42), the identity in Eq. (41) has been reduced 
to KiSAi = 0. To leading order, k = \/AM- m for each 
hole. Together with the definition of the irreducible mass 

M '» = ^ < 47 » 

where A^h is the area of the black hole's apparent hori- 
zon, we find k5A « 87r5Mi rr for each black hole. So, if 
the Komar-mass ansatz is reasonable we should find that 
J2(SMi TI )i f» along a sequence of initial-data sets. 

The second column of Table I lists the irreducible mass 
of one corotating black hole as we vary the orbital sepa- 
ration and enforce Eq. (42). Given our excision boundary 
condition, Eq. (26), the apparent horizon is, at least in- 
stantaneously, a Killing horizon. The only deviation we 
should expect in M- 1TT should be due to the fact that we 
do not have a true helical Killing vector and that the 
apparent horizon is not part of a global Killing horizon. 
As we see from the table, M; rr changes only by about 
one part in 10 4 over the entire sequence of separations. 
While small, the variations seen in M- m are larger than 
the level of accuracy of the initial data computations and 
the leading digits of the variation are also above the level 
of truncation error that result from integrating Eq. (46). 
Thus, the variation in Mj rr seems to be a true artifact of 
our approximate helical symmetry. On the other hand, it 
is quite small and this lends support to the Komar-mass 
ansatz. 



A. Effective-potential method 

Another way of testing the Komar-mass ansatz is to 
use an independent method for identifying circular or- 
bits. An "effective-potential" method for determining 
circular orbits from sequences of initial-data sets was 
outlined in Ref. [9]. This method is motivated by varia- 
tional techniques, but does not have a rigorous theoreti- 
cal foundation. Never-the-less, it has been used success- 
fully for black-hole binary initial data[9, 53, 54] and has 
been shown to agree with the Komar-mass condition in 
the case of circular orbits for thin shells of collisionless 
matter[55]. Agreement between the two methods was 
also demonstrated, within error bars, at the location of 
the inner-most stable circular orbit of sequences of circu- 
lar orbits produced using puncture data[56] (which uses 
different and more simplifying assumptions). 

The method is straightforward, but has one signifi- 
cant point of ambiguity. The effective-potential (EP) 
method identifies circular orbits as configurations with 
a minima in the binding energy along sequences of con- 
figurations where the total angular momentum is held 
fixed. The ambiguity is associated with the same free- 
dom to define the fundamental length scale mentioned 
above. In Ref. [9], this ambiguity was resolved by de- 
manding that the mass of each black hole as defined by 



TABLE I: The irreducible mass M[ TI of a single black hole 
in a corotating, non-spinning, and leading-order (LO) non- 
spinning binary as we vary the coordinate separation d. The 
length scale is set so that the ADM mass of the binary at 
very large separation is 1. Lapse boundary condition (33b) 
was used for all data. 



d 


Corotation 


M irr 

Non-spinning 


"LO" non-spinning 


40 


0.5000000 


0.5000000 


0.5000000 


35 


0.5000001 


0.5000001 


0.4999988 


30 


0.5000002 


0.5000002 


0.4999966 


25 


0.5000005 


0.5000005 


0.4999902 


20 


0.5000010 


0.5000011 


0.4999727 


19 


0.5000013 


0.5000014 


0.4999659 


18 


0.5000016 


0.5000017 


0.4999568 


17 


0.5000020 


0.5000021 


0.4999446 


16 


0.5000025 


0.5000026 


0.4999281 


15 


0.5000032 


0.5000033 


0.4999050 


14.5 


0.5000037 


0.5000037 


0.4998901 


14 


0.5000042 


0.5000042 


0.4998721 


13.5 


0.5000048 


0.5000048 


0.4998504 


13 


0.5000056 


0.5000055 


0.4998238 


12.5 


0.5000065 


0.5000064 


0.4997911 


12 


0.5000076 


0.5000074 


0.4997504 


11.5 


0.5000090 


0.5000086 


0.4996992 


11 


0.5000107 


0.5000102 


0.4996340 


10.5 


0.5000129 


0.5000120 


0.4995498 


10 


0.5000157 


0.5000144 


0.4994393 


9.5 


0.5000193 


0.5000173 


0.4992918 


9 


0.5000240 


0.5000211 


0.4990908 


8.5 


0.5000308 


0.5000259 


0.4988091 


8 


0.5000401 


0.5000330 


0.4984008 


7.5 


0.5000529 


0.5000423 


0.4977815 


7 


0.5000702 


0.5000546 


0.4967842 



the Christodoulou formula[57] be held constant along se- 
quences of configurations. 
The Christodoulou formula 

M2=M ^ r+ 4^' (48) 

includes the spin of the black hole which adds another 
level of complication. The EP method was originally ap- 
plied only to configurations in which the direction and 
magnitude of the spin was also held constant along se- 
quences of configurations. With this restriction, holding 
the Christodoulou mass M constant was equivalent to 
holding the irreducible mass M; rr constant. Furthermore, 
with this restriction on the spin, finding minima of the 
binding energy was equivalent to finding minima of the 
ADM energy. 

However, when we consider sequences of corotating bi- 
naries, the magnitude of the spins of the black holes are 
no longer held constant and we must reconsider how we 
will fix the mass freedom. From Friedman et al.[29] and 
the results shown above, is seems clear that we should 
hold the irreducible mass (or the area of the apparent 
horizon) fixed along sequences of configurations. We will 
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FIG. 1: Effective potential (EP) curves Eb/fi for corotating 
black holes vs. separation £/m. These curves are labeled by 
the orbital angular momentum J//jm which is kept constant 
along each curve. The thick red line connecting the minima of 
the EP-curves represents circular orbits; it terminates at the 
innermost stable circular orbit at the inflection point in the 
EP-curve at J/fim — 3.38. Also plotted as a dashed blue line 
is the sequence of circular orbits determined by the Komar- 
mass ansatz. 



define the total mass m and reduced mass fi as 
m = Mirrt + M irr2 , 

M = 



M irrl M irr2 



m 



and then the binding energy Eh, by 
Eh = E AnM — m 



(49) 
(50) 

(51) 



With this definition of the binding energy, minima of the 
binding energy and the ADM energy agree for corotating 
configurations as well as for sequences with the spins held 
fixed (so long as we also hold the individual irreducible 
masses fixed). 

We will, therefore, adopt the following functional defi- 
nition for the EP method. We will take as configurations 
with circular orbits, those initial data configurations that 
have a minimum of the ADM (or binding) energy along 
sequences of configurations where the total angular mo- 
mentum, irreducible masses of the black holes, and direc- 
tion of the black holes spins are held fixed. This definition 
is sufficient to handle both corotating and non-rotating 
black-hole binaries. The case of black- holes with arbi- 
trary spin will be considered in future investigations. 

Figure 1 is a plot of the dimensionlcss binding energy 
Eb/fi as a function of the dimensionless proper separa- 
tion £/m between the apparent horizons for equal-mass 
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FIG. 2: Violation of the Komar-mass condition when the ef- 
fective potential method is used to determine the sequence of 
circular orbits. Here, corotating equal-mass binaries are con- 
sidered, mfio denotes the orbital angular frequency, so that 
large separations correspond to small values of rnQo. 



corotating black holes. Each line of constant total an- 
gular momentum J/^m is an "EP curve" and the local 
minimum of that curve represents the circular-orbit con- 
figuration having that value of the total angular momen- 
tum. Passing through the set of minima, and plotted 
with a solid line, is the EP sequence of circular orbits. 
At small separations, this sequence terminates at the in- 
flection point on the first EP curve that does not contain 
a local minimum. This is the short EP curve having 
J/fim sa 3.38. The inflection point in the EP curves 
marks the termination of stable circular orbits and the 
configuration at the point is referred to as the "inner- 
most stable circular orbit" (ISCO) configuration. Also 
plotted in this figure with a dashed line is the sequence 
of circular orbits as defined by the Komar-mass ansatz. 
For a Komar sequence, the ISCO configuration is defined 
as the model with the minimum binding energy. Notice 
that the Komar and EP sequences arc nearly coincident 
except for the regime near the ISCO. 

A more quantitative comparison of the sequences ob- 
tained via the Komar-mass ansatz and the EP method 
is found by examining the error in the Komar-mass con- 
dition for circular orbits along a sequence of such orbits 
defined by the EP method. Figure 2 plots this error, 
M K — E ADM , scaled relative to the ADM energy versus 
the 5/2-power of the dimensionless orbital angular veloc- 
ity mf2o- Interestingly, we see that the relative error is 
nearly linear in (mflo) 5 ^ 2 - We see that the deviation is 
quite small, even at the ISCO near (mflo) 5 / 2 ~ 0.0032, 
where the error is roughly 0.015%. From the jaggedness 
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of the curve at large separations (small mClo), it is ev- 
ident that the errors are nearing the level of truncation 
error for the measurement of the energies. However, the 
measured deviation is clearly physical, not numerical. 

We see that the Komar-mass condition and the EP 
method appear to agree to a very high degree in deter- 
mining configurations with circular orbits. But, there is 
a measurable difference. As a final comparison, we can 
rcscale the sequence of circular orbits defined by the EP 
method so as to satisfy Eq. (42) just as was done for the 
sequence defined by the Komar-mass condition. When 
we examine the variation of M; rr along the rescaled EP 
sequence, we find the largest deviations are more than a 
factor of 100 smaller than those seen in the sequences de- 
fined by the Komar-mass condition displayed in Table I. 
This level of variability is consistent with the level of 
truncation error in the numerics. Thus, to the numerical 
precision of the calculations, we find that the sequence of 
corotating equal-mass binaries in circular orbits defined 
by the EP method satisfies the thermodynamic identity 
of Eq. (41). 



B. Spin 

Another reason for considering corotating binaries is 
that we want to calibrate our techniques for computing 
the quasi-local spin of a black hole, and corotating bi- 
naries have a physically well-defined notion of the rate 
of rotation of each black hole. From a Newtonian per- 
spective, a corotating binary has the rotational angular 
velocities of the individual objects equal to the orbital 
angular velocity. But, this is a non-local statement and, 
as such, is not well defined in the context of general rel- 
ativity. A relativistically reasonable notion of corotation 
might be to connect the rate of rotation of one black hole 
with the local rate of rotation of tidal perturbations due 
to the orbit of its companion. 

We begin by first computing the spin of each black hole 
as a function of the orbital angular velocity £Iq along a 
sequence of circular orbits that satisfy the Komar-mass 
condition. In Fig. 3, we plot the quasi-local spin as com- 
puted by Eq. (38) for the three lapse boundary conditions 
listed in Eqs. (33) and for the cases where is defined 
via the flat-space Killing vectors or via an approximate 
Killing vector as described in Sec. III. When computed 
using the flat-space Killing vectors as define by Eqs. (39), 
we note that S x and S y vanish to roundoff error. We see 
that all cases show very close agreement for the spin of 
the black holes. Also, recall that we expect our results to 
be largely independent of the choice of the lapse bound- 
ary conditions and plot the results from different lapse 
boundary conditions to confirm this conjecture. 

While the various measures of the spin agree well, we 
have yet to determine if the magnitude of the spin is 
correct. To do so, we first note that the spin of a Kerr 
black-hole, Skcm is given in terms of its irreducible mass 
Mj rr and the rotational angular frequency of its horizon, 
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FIG. 3: Quasi-local spin S/M^ rT of a black hole in a corotating 
equal-mass (CO) binary along the sequence of circular orbits 
(parameterized by the orbital angular velocity milo). S z and 
Sk are denned in Eq. (40). Results are given for sequences 
with three different lapse boundary conditions. 



Sl B , by 



SKerr(-Mirr) 



sjl - 4(A/ irr fb) 2 



(52) 



So, if we know the rotational angular frequency of the 
black hole's horizon, then we can check how well S z and 
Sk satisfy Eq. (52). Conversely, by inverting Eq. (52) we 
can estimate the angular velocity Sl B of the black hole in 
our numerical initial data from the measured quasi-local 
spins. It will be convenient to express this as 



r> 



i 



S/M'i 



Mi rr mQo 



s]l + \{S/Ml) 



(53) 



In Ref. [5], we assumed that the rotational angular 
velocity of the black holes in a corotating binary was 
equal to the orbital angular velocity fl B = ^o- We will 
refer to this assumption as the "leading-order" estimate. 
The lower six lines of Fig. 4 display the relative error 
between the measured quasi-local spin and the spin of a 
Kerr black hole, Eq. (52), with the same irreducible mass 
and with £l B = ^o- We see that the quasi- local spin 
estimate derived from the approximate Killing vectors is 
slightly better than the estimates based on the flat-space 
Killing vectors. Furthermore, the variation due to the 
lapse boundary condition is smaller for the approximate 
Killing vector cases. In all cases, the relative error is 
smaller than about 8% up to the ISCO which occurs near 
mflo s=» 0.1, a quite large disagreement. 

In order to make a better estimate of the error in the 
quasi-local spin measurements, we need a better estimate 
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FIG. 4: Difference between the quasi-local spin of a black 
hole in a corotating (CO) binary and the result of the Kerr- 
formula Eq. (52). For the lower set of six lines (Qb = flo), 
the rotation rate of the black hole is simply taken to be equal 
to the orbital angular velocity. For the upper set of six lines 
(Qb = Qt), the rotation rate is taken to be equal to that 
of the tidal field of the companion hole as measured in the 
LARF. Each set of six lines corresponds to the cases plotted 
in Fig. 3. 
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FIG. 5: Ratio of the rotational angular velocity SIb of a black 
hole, as determined from the quasi-local spin Sk, to the or- 
bital angular velocity of the binary along the corotating (CO) 
sequence. The dotted line (labeled IPN) represents the lead- 
ing order correction to the analytic estimate of the ratio. The 
dashed line (labeled 1PN+) shows the higher order correction 
of Eq. (56) with A = T = 0. Numerical results are given for 
three different lapse boundary conditions. 



of the correct spin value. To improve on this estimate 
we note that Alvi[58] has computed the leading order 
correction to the rotation rate Ot of the tidal field of a 
companion star as measured in the local asymptotic rest 
frame (LARF) [59] of a black hole 



Qj 1 — £1q 



(54) 



where rj = fi/m and b is the separation of the black holes 
in harmonic coordinates. We can express this in a gauge 
independent way using a post-Newtonian expansion for 
b/m obtained for circular orbits [60] 



? = («^o) 2/3 





1 + (1 - ^f]){mn ) 2/3 + o{mn ) 



Substituting, we find 



(55) 



On 



= 1 -r](mn ) 2/3 + A(mOo) 



(56) 



[ v (i - |r?) - r] (mn ) 4 / 3 + O(mO )t, 



where A and T are functions of rj coming from the un- 
known terms of order (m/6) 3 / 2 and (to/6) 2 respectively 
in Eq. (54). The leading order error term in Eq. (55), 
which includes spin-orbit coupling terms, contributes to 
the term of order (toOq) 2 in Eq. (56). We now make 



the physical assumption that fix represents the angu- 
lar velocity fls at which a black hole should rotate in 
its LARF in order to be in corotation. This assumption 
improves the agreement with the Kerr-formula dramati- 
cally, as shown by the upper set of six lines in Fig. 4. 

Figure 5 shows, from a different perspective, the 
improvement obtained by using Alvi's result. From 
Eq. (53), we obtain the expected value of f2s/fin, from the 
initial data, assuming the Kerr-formula is exactly satis- 
fied. In Fig. 5, we show this expected value for flg/fl for 
each of the three lapse boundary conditions and compare 
these to fir/fin. from Eq. (56). The dotted line labeled 
IPN includes only the leading-order correction to CIt/Ho- 
Since we are plotting the ratio against (toOo) 2 ^ 3 , this 
line is linear. The dashed line labeled 1PN+ shows the 
full analytic estimate assuming the unknown coefficients 
A = r = 0. For our equal-mass binaries, 77 = 1/4. 

We see that the analytic estimate for the rotation rate 
is in reasonably good agreement with the numerical re- 
sults. In fact, the difference between them is nearly lin- 
ear in mfin. and is well fit by setting A ~ — 0.085(for the 
equal-mass case of 77 = 1/4 and T = 0). It would be quite 
interesting to determine the next term in the rate of ro- 
tation of the tidal field given in Eq. (54) since this term 
would fix A. We should note, however, that we cannot 
expect our estimates of the spin to agree perfectly with 
the higher-order analytic estimates of the black-hole ro- 
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tation rate for corotation. There are two reasons for this. 
One is that our initial data is conformally fiat. In Ref.[5], 
we examined the case of a single rotating, conformally- 
flat black hole. Maximum rotation in a corotating binary 
occurs at ISCO where the rotational angular velocity has 
a value of M 1ti {Ib ~ 0.05. At this rate of spin, we expect 
the approximation of conformal flatness to introduce an 
error of about 0.02%. The second source of error is the 
quasi-local estimate of the spin itself which has some a 
priori unknown level of uncertainty. We show a rough 
indication of the current level of uncertainty in Fig. 4. In 
the top six lines, we again plot the relative error of the 
quasi-local spins, but now plotting the error relative to 
the Kerr spin based on the tidal rotation rate given by 
Eq. (56) with A = T = 

We see that at ISCO (mfio <~ 0.1), the relative errors 
for the spins computed using an approximate Killing vec- 
tor are all less than 1%, but significantly larger than the 
0.02% error caused by the assumption of conformal flat- 
ness. Unfortunately, we cannot determine whether this 
remaining error is due to the inherent uncertainty of our 
quasi-local measure of the spin or to the unknown value 
of A (and higher-order terms). The only thing that we 
can conclude is that the uncertainty in the spin based on 
comparing the two different approaches to compute the 
spin is of the order of 1%, and is about 0.5% based on 
the uncertainty introduced by the different lapse bound- 
ary conditions. Both are of the same order as the total 
error which is remarkably small and gives us considerable 
confidence in our quasi-local spin measure. 

V. NON-SPINNING BINARIES 

In order to model non-spinning binaries, we must 
choose an appropriate non-vanishing excision boundary 
condition on /3y\ Since our models are set up so that the 
orbital angular momentum is pointed in the positive z 
direction, we choose the boundary condition on (3^ via 
Eq. (32) and set = Q z y the flat-space Killing vector of 

Eq. (39c). This choice for £ l is appropriate because we 
have a spherical excision surface in a flat conformal geom- 
etry, so is a conformal Killing vector of the physical 
induced metric on the excision surface. The only freedom 
remaining in the boundary condition is the magnitude of 

A. Leading-order approximation: fl r = fio 

In Ref. [5], it was argued that, at least to leading or- 
der, the correct choice for non-spinning black holes was 
to choose £l r to be equal to the orbital angular velocity, 
f2 r = Qq. The numerical results were in good agreement 
with post-Newtonian estimates for sequences of equal- 
mass non-spinning binaries in circular orbits and with 
post-Newtonian estimates for the location of the ISCO. 
In this section, we revisit this argument with the new 
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FIG. 6: EP curves Eb/fi for equal-mass "leading-order" non- 
spinning binaries (fi r = fio) plotted vs. separation £/m. These 
curves are labeled by the value of J//xm along each curve. 
Also plotted are the line connecting the minima of the EP- 
curves, as well as the sequence of circular orbits as determined 
by the Komar-mass ansatz. The Komar-mass ansatz and the 
effective potential method clearly disagree. 

diagnostic tools presented earlier in this paper, in partic- 
ular, effective potential plots and quasi-local measures of 
the spin. 

Figure 6 displays several EP curves for the leading- 
order (i.e. fi r = fio) non-spinning binaries where the 
total angular momentum is held fixed. It is immediately 
clear that the sequence of circular orbits determined by 
the Komar-mass condition does not intersect the minima 
of the EP curves. Furthermore, there does not appear 
to be an inflection point in the EP curves anywhere near 
the minimum of the Komar sequence of circular orbits. 
This clearly indicates that a problem of some kind exists. 

The source of the problem is made clear by examining 
the quasi-local spins of black holes in our leading-order 
non-spinning models. The inset of Fig. 7 shows the mag- 
nitude of the spin on one of the black holes in the leading- 
order non-spinning equal-mass binary as a function of the 
orbital angular velocity. As in the corotating case, the x 
and y components of the spin vanish to roundoff error and 
we see that the z component of the spin computed using 
flat-space Killing vectors agrees well with the magnitude 
of the spin computed using an approximate Killing vec- 
tor. But, it is clear that the magnitude of the spin is not 
zero. To get a better understanding of whether or not 
this magnitude of spin is small, Fig. 7 also compares this 
spin with the spin of a Kerr black hole having a rotational 
angular velocity equal to the binary's orbital angular ve- 
locity. Essentially, this is a measure of the magnitude of 
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FIG. 7: Quasi-local spin of the "leading order" non-spinning 
(LN) binaries (as shown in Fig. 6) plotted vs. orbital angu- 
lar velocity. The inset displays the dimensionless spin value 
S/M? r . The main figure displays the magnitude of the spin 
relative to the spin of a Kerr black hole with Qb = f2o- The 
residual spin is quite large. 



the spin relative to the corresponding corotating case. 

We see that the spin for the leading-order non-spinning 
equal-mass binaries can be as large as 30% of the corota- 
tion spin at ISCO. This is clearly a significant deviation 
from the desired non-spinning configuration we wish to 
model. 

Another indication that there is an error in the leading- 
order non-spinning equal-mass binaries is found in exam- 
ining how well the sequence of leading-order non-spinning 
circular orbits agrees with the thermodynamic identity of 
Eq. (41). In Table I, we show the results for sequences of 
circular orbits constructed using the Komar-mass condi- 
tion. The last column shows the results for leading-order 
non-spinning binaries. We find the disturbing result that 
the irreducible mass decreases as the binary evolves to 
smaller radii. More importantly, comparing this varia- 
tion to that seen in the second column, we find that the 
size of the variation in M; rr for the leading-order non- 
spinning case is roughly 20 times larger than that seen 
in the corotating case (see also Fig. 9 of Ref. [5]). Of 
course, the thermodynamic identity was derived for coro- 
tating binaries, not non-spinning binaries so we shouldn't 
be surprised that it is not satisfied in this case, and by 
itself, this failure is not too worrisome. However, the 
combined evidence in this section, in particular Figures 
6 and 7, makes it clear that the leading order approxi- 
mation Cl r — Oo is not satisfactory. Since we can now 
evaluate the quasi-local spin, we no longer have to choose 
£l r by some ad hoc prescription - instead, we can simply 
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1/m 

FIG. 8: EP curves Eb/n for "true" non-spinning black holes 
vs. separation t/m. The curves are labeled by the orbital 
angular momentum J/fj,m which is kept constant along each 
curve. The thick red line connecting the minima of the EP- 
curves represents circular orbits, and terminates at the in- 
nermost stable circular orbit at the inflection point in the 
EP-curve at J/fim — 3.12. Also plotted as a dashed blue 
line is the sequence of circular orbits as determined by the 
Komar-mass ansatz. 

choose it such that the quasi-local spin vanishes. This 
approach will be explored in the following section. 

B. Correct approach for non-spinning binaries 

To correctly model non-spinning black-hole binaries, 
we need to choose Q r so that the quasi-local measure of 
the spin vanishes, Sk = 0. This will involve root-finding, 
i.e., the constraints will have to be solved for different 
values Q r until the solution satisfies Sk = 0. Since the 
correct value of il r for non-spinning black holes will be 
close to fio, we define a rotation fraction f r as 

Q r = f r Q - (57) 

While the definition of f r mirrors Eq. (56) from our com- 
parison with post-Newtonian results, we note that those 
post-Newtonian results do not enter the construction of 
non-spinning binary black hole initial data. In practice, 
we determine f r so that the quasi-local spin based on the 
approximate Killing vector vanishes to about one part in 
10 8 or better. 

In Fig. 8 we first plot the EP curves for "true" non- 
spinning equal-mass black-hole binaries. The solid line 
passing through the minima of the EP curves is the se- 
quence of circular orbits as defined by the EP method. 
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FIG. 9: Quasi-local spin for a "true" non-spinning (NS) black 
hole in an equal-mass binary along the sequence of circular 
orbits, parametrized by the orbital angular velocity. The same 
six cases are plotted as in the inset of Fig. (7). Note that here 
Sk = by construction. 
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FIG. 10: Rotation parameter f r = Q r /Q,o along the sequence 
of non-spinning (NS) black holes, parametrized by the 2/3- 
power of the orbital angular velocity. The dotted line plots 
the fit to f r through the term of order (mQo) 2 ^ 3 . 



This sequence terminates at small £/m at the ISCO which 
is defined by the occurrence of an inflection point in an 
EP curve. Finally, plotted as a dashed line in this fig- 
ure is the sequence of circular orbits as defined by the 
Komar-mass condition. It is immediately apparent that 
the EP curves have the correct qualitative behavior and 
that the Komar sequence of circular orbits, drawn as a 
dashed line, passes very close to the minima of the EP 
curves everywhere except near the ISCO. 

In Fig. 9, we plot the value of the quasi-local spin of one 
black hole in an equal-mass binary as a function of the 
orbital angular velocity. The sequences are circular or- 
bits defined by the Komar-mass condition and we include 
plots of the three different lapse boundary conditions for 
both definitions of the spin. By construction, the quasi- 
local spin defined by the approximate Killing vector is 
zero, Sk = 0. However, the quasi-local spin defined by 
the flat-space Killing vector (S z ) is not necessarily zero 
and is a rough measure of the uncertainty in our defini- 
tion of the spin. For separations larger than ISCO, we 
see that, in terms of S z , the spin does not exceed 0.4% 
of maximal rotation. 

In Fig. 10, we plot the value of the rotation fraction 
f r = Q r /Qo as a function of (mfio) 2 ^ 3 - To leading order, 
we expected f r = l and that to model non-rotating black 
holes, we needed to set the excision boundary conditions 
on 1 so as to "unrotate" the black hole from the case of 
corotation when (3^ = 0. We might then expect that at 
higher order we would find f r = fls/^o ~ ^t/^o- How- 
ever, comparing f r in Fig. 10 to fis/fio in Fig. 5 we see 



that f r is significantly smaller. In terms of the previous 
physical interpretation of the excision boundary condi- 
tions for non-spinning holes, this implies that we need 
to "unrotate" the black hole by less than the corotation 
rate. Currently, we do not have a theoretical interpreta- 
tion for the value of f r seen in Fig. 10. If we assume that 
an expansion in terms of mf2o takes the same functional 
form as flx/^o m Eq. (56), we find that f r is well fit by 

/,. = l-0.925(mfi ) 2/3 +0.36(mr2o)-l-4(mr2 ) 4/3 - (58) 

We note, however, that f r can also be well fit by a func- 
tion that includes a term of order (rnflo) 1 ^ 3 - 

C. Results for non-spinning binaries 

In Ref. [5], we carefully examined the case of equal- 
mass non-spinning binaries, but the initial-data sets were 
constructed using the method for defining non-spinning 
black holes that is correct only to leading order. Be- 
cause our improved approach does yield a different spin 
for the resulting black holes, we reexamine the physical 
content of these configurations. In the previous section, 
we computed the sequence of circular orbits for truly non- 
spinning binaries (cf. Fig. 8). In Figs. 11-13, we plot 
parameters along the sequence of non-spinning binaries 
defined using the Komar-mass condition. In each case, 
we compare our initial-data sequences to the effective one 
body (EOB) post-Newtonian results[5, 10] at first, sec- 
ond, and third post-Newtonian order, to the results of an 
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FIG. 11: Binding energy vs. total angular momentum along 
the sequence of non-spinning (NS) equal-mass black holes 
(as defined by the Komar-mass condition). For comparison, 
the leading-order non-spinning (LN) results from our earlier 
work [5] are included, as well as results from Refs. [9, 10]. 

earlier initial-data method [9], and to the results based on 
the leading-order method of defining non-spinning black 
holes. While the difference between our improved numer- 
ical results and those based on the leading-order method 
are not dramatic, they do become significant at small 
separations. 

The most remarkable change produced by the im- 
proved method is seen in Fig. 11 which plots the dimen- 
sionlcss binding energy as a function of the dimension- 
less total angular momentum. Using the leading-order 
method, the sequence did not approximate a cusp at the 
ISCO as we would have expected. However the improved 
data clearly approximates a cusp. In Ref.[5], we pointed 
out that we did not understand why our non-spinning 
data lacked this feature. We now understand that the 
approximate cusp is a necessary feature of a sequence 
that is in good agreement with the EP method. This can 
be seen by looking at Figs. 6 and 8 and considering the 
behavior near the inflection point of Fig. 8 that defines 
the ISCO for the EP method. In fact, a sequence of cir- 
cular orbits defined by the EP method necessarily has 
the minima in Eb/ji and J/fim coincide, resulting in an 
exact cusp in the sequence at ISCO. 

We can also reexamine how well the improved non- 
spinning data agrees with the thermodynamic identity 
of Eq. (41). The third column of Table I shows M irr 
along a sequence of non-spinning equal-mass binaries in 
circular orbit constructed using the Komar-mass condi- 
tion. Comparing to the results of the second column for 
the corresponding case of corotation, we see that the new 
approach for defining non-spinning binaries yields results 



FIG. 12: Binding energy vs. orbital frequency along the se- 
quence of non-spinning (NS) equal-mass black holes (as de- 
fined by the Komar-mass condition). Lines are labeled as in 
Fig. 11. 
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FIG. 13: Total angular momentum vs. orbital frequency 
along the sequence of non-spinning (NS) equal-mass black 
holes (as defined by the Komar-mass condition). Lines are 
labeled as in Fig. 11. 



that are comparable in magnitude. Also, the variation in 
Mj rr is much smaller than those seen in leading-order 
non-spinning data of the last column. 

Finally, it is interesting to reconsider the sequence 
of non-spinning circular orbits constructed by the EP 
method. Recalling that an exact stationary solution of 



16 




FIG. 14: Violation of the Komar-mass condition when the 
effective potential method is used to determine the sequence 
of circular orbits. Here, non- spinning (NS) equal-mass bina- 
ries are considered, and the sequence is parametrized by the 
orbital angular momentum. 



Einstein's equations for a black-hole binary can only be 
found if the black holes are in corotation, we should ex- 
amine again how well the Komar-mass condition and the 
EP method agree in their predictions of circular orbits. 
From Fig. 8, we see that the sequences of circular or- 
bits denned by the two methods nearly coincide except 
for the regime near the ISCO. For a more quantitative 
comparison, we can again examine the relative error in 
the Komar-mass condition. Figure 14 displays this error 
for the case of non-spinning equal-mass binaries plotted 
against the dimcnsionlcss orbital angular velocity. The 
magnitude of the relative error is comparable to that of 
the corotating case plotted in Fig. 2. 

We have also examined how well the sequence of non- 
spinning circular orbits constructed by the EP method 
agrees with the thermodynamic identity of Eq. (41). As 
with the corotating case, we use the freedom to define the 
fundamental length scale along the sequence to enforce 
Eq. (42). Then the deviation in M- lTI along the sequence 
is a measure of how well the thermodynamic identity is 
satisfied. As with the case of corotation, we find that 
the variation in the mass is smaller along the EP se- 
quence than in the Komar-mass sequence. However, the 
variations are not as small as seen in the corotating EP 
sequence. While the variation in the corotating EP se- 
quence were consistent with truncation error, the varia- 
tions in the non-spinning EP sequence are consistently 
larger by a factor of about five and the increase in the 
mass appears to be systematic and significant. 

This result should not be surprising since the ther- 



modynamic identity (41) we are testing has only been 
defined for the case of corotating binaries. What is re- 
markable is that this identity is satisfied so well for the 
case of non-spinning black holes. 

VI. DISCUSSION 

Our purpose in this paper has been to explore the spins 
of black holes in equal-mass binaries in order to verify 
that the corotating and non-spinning cases were being 
modeled correctly, and also to explore the assumptions 
being used to identify configurations that are in circu- 
lar orbits. In the process of doing this, we have dis- 
covered that the assumptions made in the first attempt 
to construct non-spinning equal-mass binaries (using the 
quasicquilibrium approach described in this paper) [5] 
were not leading to sufficiently accurate representations 
of non-spinning binaries. However, the same quasi-local 
techniques used to measure the spins of the black holes 
can also be used to define a new approach for construct- 
ing non-spinning black-hole binaries, and this approach 
has produced excellent results. A detailed description of 
the new results is found in Sec. VC and in App. C. 

TABLE II: Parameters of the ISCO configuration for non- 
spinning equal-mass black holes. Results for the ISCO are 
given for three different choices of the lapse boundary condi- 
tion with circular orbits defined by the Komar-mass ansatz 
and for a single lapse boundary condition with circular orbits 
defined by the EP method. For comparison, the lower part of 
the table lists results of Refs. [9-11]; "Conf. Imag." [9] repre- 
sents data derived from previous numerical initial-data sets; 
"PN standard" [11] represents a post-Newtonian expansion in 
the standard form without use of the EOB-technique. 



Lapse BC ISCO type 


m Qq 


E b /m 


J/m 2 


iW) _ o Komar 


0.122 


-0.0194 


0.779 


d(^) = cap Komar 

clr Zr 

EP 


0.122 
0.121 


-0.0194 
-0.0193 


0.779 
0.780 


aip = 5 Komar 


0.124 


-0.0195 


0.778 


Conf. Imag. 


0.166 


-0.0225 


0.744 


IPN EOB 
2PN EOB 
3PN EOB 


0.0692 
0.0732 
0.0882 


-0.0144 
-0.0150 
-0.0167 


0.866 
0.852 
0.820 


IPN standard 
2PN standard 
3PN standard 


0.5224 
0.1371 
0.1287 


-0.0405 
-0.0199 
-0.0193 


0.621 
0.779 
0.786 



For completeness, we also include tables and figures 
that detail the ISCO configurations for both corotating 
and non-spinning cases. Table II displays the new results 
for ISCO configurations for non-spinning equal-mass bi- 
naries. These have changed considerably as is evident 
by examining these same parameters plotted in Figs. 15— 
17. In these figures, the previous non-spinning results 
are labeled as LN : QE and the new results as NS : QE. 
To facilitate direct comparison, we also include Table III 
which displays the ISCO result for corotating equal-mass 
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binaries. In both cases, wc include results obtained using 
the three choices for the lapse boundary condition given 
in Eq. (33) and with circular orbits determined using the 
Komar-mass condition. Both tables also include a single 
ISCO model where the circular orbit is determined us- 
ing the EP method. Wc include the results only for the 
lapse boundary condition in Eq. (33b) because locating 
an ISCO model within the EP approach is very expensive 
computationally and the results are not sufficiently dif- 
ferent from those obtained via the Komar-mass condition 
to warrant the expense. 



TABLE III: Parameters of the ISCO configuration for corotat- 
ing equal-mass black holes. Results for the ISCO are given for 
three different choices of the lapse boundary condition with 
circular orbits defined by the Komar-mass ansatz and for a 
single lapse boundary condition with circular orbits defined 
by the EP method. Layout as in Table II. 



Lapse BC ISCO type 




E b /m 


J/m 2 


<H««j _ Komar 


0.107 


-0.0165 


0.844 


d(^) _ o* Komar 
EP 


0.107 
0.104 


-0.0165 
-0.0163 


0.843 
0.845 


aip = i Komar 


0.107 


-0.0165 


0.843 


HKV-GGB 


0.103 


-0.017 


0.839 


IPN EOB 
2PN EOB 
3PN EOB 


0.0667 
0.0715 
0.0979 


-0.0133 
-0.0138 
-0.0157 


0.907 
0.893 
0.860 


IPN standard 
2PN standard 
3PN standard 


0.5224 
0.0809 
0.0915 


-0.0405 
-0.0145 
-0.0153 


0.621 
0.882 
0.867 
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FIG. 15: ISCO configuration for non-spinning (NS) binary 
black holes, computed with three different lapse boundary 
conditions. For comparison, the leading-order non-spinning 
(LN) and corotating (CO) binary black holes from our earlier 
work [5] are included, as well as results of Refs. [8-11]. For 
post-Newtonian calculations the size of the symbol indicates 
the order, the largest symbol being 3PN. "PN standard" [11] 
represents a PN-expansion in the standard form without use 
of the EOB-technique (only 2PN and 3PN are plotted). 



We think it worth noting that the value in consider- 
ing the ISCO parameters is questionable. The ISCO, of 
course, can only be defined when the effects of radiation 
reaction are ignored so that the dynamics is conservative. 
The effects of radiation reaction are certainly not small 
for equal-mass binaries near ISCO, and so the physical 
value of locating the ISCO is limited at best. However, 
the ISCO is a well-defined, unique point in a sequence 
of (quasi-) circular orbits and so has value as a point of 
comparison between various methods. It also has value 
in that it marks the regime where the quasiequilibrium 
approximations have clearly broken down. 

In addition to yielding improved initial data for non- 
spinning equal-mass black-hole binaries, our investiga- 
tions have shown that the application of quasi-local 
methods for determining the spins of black holes in bi- 
nary systems yield remarkably good results. For the case 
of corotating binaries, we have shown that the quasi-local 
spin measured for the individual black holes is in excellent 
agreement with the theoretical expectation for the spin. 
However, a nagging question remains with regard to the 
computation of quasi-local spins. This question is asso- 
ciated with the meaning of "approximate solutions" of 
Killing's equation when no true symmetry exists. This is 
an issue we hope to examine further. We note that during 
the preparation of this manuscript, Schnettcr et al.[61] 
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FIG. 16: ISCO configuration for non-spinning (NS) binary 
black holes, computed with three different choices of the lapse 
boundary condition. Symbols as in Fig. 15. 
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FIG. 17: ISCO configuration for non-spinning (NS) binary 
black holes, computed with three different choices of the lapse 
boundary condition. Symbols as in Fig. 15. 



Perhaps more important than the issue of which ap- 
proach is better is the fact that both methods agree so 
well also for the case of non-spinning binaries and that 
both methods are in good agreement with the thermody- 
namic identity of Eq. (41). This is surprising because the 
identity was intended to apply only in the case of coro- 
tating black- holes binaries. If the agreement only held at 
large separation, it would be reasonable to assume that 
the identity held because strong-field effects were simply 
small in this regime. However, the agreement holds at all 
separations which suggests that there may be something 
more fundamental at work. 
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have also raised this issue. 

Concerning the theoretical expectation for the spin, we 
note that the leading order result, obtained by assuming 
that a corotating black hole rotates with an angular ve- 
locity equal to the orbital angular velocity of the binary, 
is not adequate. Using Alvi's[58] leading-order correction 
to the rate at which a tidal distortion produced by an or- 
biting star travels around a black hole, we have deduced 
an improved prediction for the spin of a corotating black 
hole. It would be useful if a higher order calculation of 
Alvi's tidal distortion rotation rate were available. The 
effect of the corrected rotation rate on post-Newtonian 
models of corotating binaries should also be considered, 
especially when higher order spin-orbit and spin-spin in- 
teraction terms become available. 

Finally, we note that there is remarkably good agree- 
ment between sequences of circular orbits defined via the 
Komar-mass ansatz and by the EP method. It is natural 
to ask which method should be used to locate configu- 
rations in circular orbit. As described in Sec. IV A, we 
find that the thermodynamic identity (41) of Friedman 
et al.[29] that applies to corotating binaries in circular 
orbits is satisfied much more closely by models defined 
by the EP method than those that satisfy the Komar- 
mass ansatz. While the EP method yields deviations 
in the one free quantity {M- m ) that are two orders of 
magnitude smaller than the Komar-mass condition, the 
deviations produced by either approach are very small. 
From a practical perspective, applying the EP method 
is quite expensive computationally since it requires find- 
ing the minima of a "function" that is itself defined by 
root-finding methods. For most applications, it is hard 
to justify this additional cost. 



APPENDIX A: SOLVING KILLING'S EQUATION 

A straightforward method for finding Killing vectors 
was outlined by Dreyer et al.[27]. Here we outline that 
method as we have implemented it to locate Killing vec- 
tors on a surface S with S 2 topology. In our case, S 
is a coordinate sphere in the flat conformal 3-geomctry. 
Following Eq. (27), our metric is denoted hij and it is 
conformally related to which is the metric of a coor- 
dinate 2-sphere with radius r. Following Ref. [27], we let 
£ l denote a Killing vector on (S,hij) and define a two- 
form Lij = Di£j, where Z), is the covariant derivative 
compatible with hij. Now, the Killing equation can be 
written simply as 

L {ij) = 0. (Al) 

However, instead of solving the Killing equation directly, 
Dreyer et al. propose to solve the Killing transport equa- 
tions: 

D£j = L^, (A2a) 
DiLjk = 2 R k ji%, (A2b) 

for the two-form and vector £\ Here, 2 Rkji l is the 
Ricmann tensor associated with hij . Note that Eq. (A2b) 
follows from Eq. (A2a) by using various symmetries of 
Riemann and assuming that £* is a Killing vector (i.e. 
that D(i£j) = 0). 

The central idea in using Eqs. (A2) to find Killing vec- 
tors is to note that the Killing transport equations act 
as a linear map on the set of variables Lij). In par- 
ticular, consider integrating these equations along some 
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path from point p to point q. If {S,h%j) has solutions 
of the Killing equation (Al), and if Lij)\ p come from 
one of these Killing solutions, then integrating the Killing 
transport equations along the path will yield (£,i,Lij)\ q 
which comes from the same Killing vector. 

Since our surface S is 2-dimensional, it follows that 



(A3) 
(A4) 



where 2 R is the 2-dimensional Ricci scalar associated with 
hij and L is an oriented scalar on <S. For our particular 
metric, 



ds 2 = i/)V(dr + sin 2 e&cj) 1 ), 



(A5) 



the Killing transport equations can be written in spheri- 
cal coordinate components as 



^ = 2£ e d e blip - -r^-^fy. In?/', 
30 sm z e 



(A6a) 



56> 



3L 



cot 0^ + Lip 4 r 2 sin + 2£ e d 4 , In tJj + 2^3 e In ip, 

(A6b) 

■IirD/.v.r) -^i'... (A6c) 



sin 



cot0^ - Lip 4 r 2 sin6» + 2£ e <9 lnV> + 2^d e Imp, 

(A7a) 



^j- — - sin 6> cos - 2 sin 2 0^gde \nip + 2^3^ In ?/>, 



3L ( 1 



2D l D l sm9i 



(A7b) 
(A7c) 



It is particularly convenient to compute these quanti- 
ties in terms of a basis set (0, (f) which is orthonormal as 
defined on the unit 2-sphere. We find then that: 



= & and 6 



1 



sin 



(A8) 



The Killing transport equations then take the form 



30 



= 2&(S71nV0s-2&(571nV0j 



(A9a) 



-j± = LV/r 2 + 2£ e -(FlnV^ + 2^( s Vln^ e -, (A9b) 
_ (l-2y 2 ln^) 
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(A9c) 



cos 6>^ - Lifj 4 r 2 sin 6» (AlOa) 

+ 2 sin %( S V In ^ + 2 sin ^(J7 In <%, 

- cos % - 2 sin %(y In + 2 sin 0^(.V In 

(AlOb) 



0& 



<9L _ (l-2 b V 2 hn/>) , 



t\) A r 2 



sin c?£g 



(AlOc) 




where y and J7 2 are the usual gradient and Laplacian 
operators defined for the unit 2-sphcrc. 
Let V denote the vector of quantities 



(All) 



Using fourth-order Runge-Kutta, we integrate Eqs. (A10) 
around a closed path from <p ~ to <j) = 2tt along the 
equator = -| starting with three different values for 
V| w=0) : (1,0,0), (0,1,0), and (0,0,1). The three result- 
ing vectors, V|(0 =2 tt) can be used to construct a matrix 
M that represents the action of the linear map of the 
Killing transport equations on any vector: 



VI 



M- V 



(A12) 



If V|(^, =0 ) is derived from a Killing vector, then Eq. (A12) 
will yield V|(0 =2 7r) = V| ( ^ =0) . So, if (S,hij) possesses a 
global solution of the Killing equation, it will be asso- 
ciated with a unit eigenvalue of M and the associated 
eigenvector will be derived from that Killing vector. In 
general, (iS, hij) will not possess a Killing vector and M 
will have no unit eigenvalues. One could proceed to con- 
struct an "approximate Killing vector" by using an eigen- 
value that is sufficiently close to unity. 

However, for the corotating and non-spinning black- 
hole binary initial data we consider in this paper, the 
conformal factor possesses a reflection symmetry through 
the plane of the orbit. We thus find that the eigenvec- 
tor (0, 1, 0) always has unit eigenvalue when we integrate 
along the equator. 2 Unfortunately, this vector does not 
necessarily represent a global Killing vector on S. 

Our pseudo-spectral code for solving the constraint 
equations represents quantities on the excision bound- 
aries in terms of spherical harmonic decompositions. To 
proceed we need to populate the collocation points (grid 
points) of the spectral grid with values for V obtained by 
propagating the eigenvector (0, 1, 0) at (0 = 5-, (j) = 0) to 



2 This is most easily seen by writing the Killing transport equa- 
tions in terms of an orthonormal basis set (6, </>) denned with 



respect to the full metric hij so that 



ii'-'r sin 



5e and § r 
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all of these points via the Killing transport equations. If 
our starting vector were constructed from a global Killing 
vector, then the path we take to each of these points 
would not matter. However, since we do not in general 
have a true global Killing vector, we must specify how 
we populate the grid points. 

First, we integrate Eqs. (A10) to populate values along 
the equator (# = §) wherever they are needed, integrat- 
ing from 4> = to 4> = 27r. For each line of colatitude con- 
taining collocation points, we integrate Eqs. (A9a-A9c) 
starting with the known values at the equator. All inte- 
grations are performed using fourth-order Rungc-Kutta. 
With values for ^ and £i at all collocation points we can 
construct a vector spherical harmonic representation of 
our solution over the entire surface S. All that remains 
is to normalize the Killing vector and check to see if we 
have a true Killing vector. 

A rotational Killing vector (as we are trying to con- 
struct) is normalized so that its affine length is 27r. Con- 
sider a path parameterized by t and defined by 



1 



m 



dt ^ 4 r 2 



1 



dt ip 4 r 2 sm8 



(Al3a) 
(A13b) 



We integrate from some starting location (0,4>)\o with 
t = until the path closes at the starting point. If £ is 
not properly normalized, then the final value of tf will not 
be 2ir. However, using tf we can rescale £ so that it does 
have an affine length of 2-7T. In practice, we normalize the 
solution by integrating along the equator, then check that 
the solution is correctly normalized to within truncation 
error for several different integral paths starting at (f> = 
and various initial values for 9. 

There are many possible ways of determining whether 
or not our solution is a global solution of the Killing 
equation. For simplicity, we check the following set of 
scalar conditions that must be satisfied everywhere if our 
solution is a true Killing vector of (5, hij): 

CD t L = 0, (Al4a) 
AC = 0, (A14b) 
DX, = 2L. (A14c) 

In terms of our basis defined by Eq. (A8), these become 



ip 4 r 



|-a (ts$L)s + Z$(pL)j) = 0, (Al5a) 



1 



(A15b) 



-^px£=2L, (A15c) 



and y- and J/x are the usual divergence and twist op- 
erators defined for the unit 2-sphere. We find that these 
identities arc never satisfied for our solutions although 



the residuals seem "small" and decrease in size as the 
separation between the black holes increases. Because 
we have not yet found any way to meaningfully normal- 
ize the values of the residuals, we do no bother to report 
their values here. 

We do note that, because we represent £ in terms of 
vector spherical harmonics, it is possible to filter the co- 
efficients in the vector spherical harmonic expansion to 
ensure that it is divergenceless. That is, wc can, a pos- 
teriori guarantee that Eq. (Al5b) is satisfied. However, 
this does not change the other diagnostics. All of the 
results presented in this paper that depend on solutions 
of the Killing equation have been filtered in this way. 
We note that this filtering has a negligible effect on all 
measured quantities. 

We conclude this Appendix by noting that our solu- 
tions of the Killing equation are "approximate Killing 
vectors" in some sense that is not at all well defined. In 
all cases, we find that our approximate Killing vector is 
very similar to the corresponding solution on the unit 2- 
spherc. These are conformal Killing vectors of (S,hij) 
as written down in Eqs. (39a-39c). We leave it to fu- 
ture work to more rigorously define the meaning of these 
"approximate Killing vectors." 

Finally, we note that during the preparation of this 
manuscript, Schnetter et al.[61] have shown that when 
a divergence free approximate Killing vector is used to 
define the quasi- local angular momentum, as in Eq. (38), 
then it is gauge invariant. Therefore, they suggest that 
such a divergence-free approximate Killing vector "can 
be viewed as an ersatz axial symmetry vector even in the 
absence of axisymmetry." 



APPENDIX B: COROTATING SEQUENCE 

In this Appendix, we list the numerical results for coro- 
tating equal-mass black holes in a quasi-circular orbit as 
defined by the Komar ansatz. We have assumed confor- 
mal flatness, maximal slicing, and used Eq. (33b) for the 
lapse boundary condition on both excision surfaces. The 
data is scaled relative to the total and reduced masses 
(m and /x) defined with respect to the irreducible mass 
of the apparent horizons M lvr . 

In Table IV, d and Qq are, respectively, the separa- 
tion of the centers of the excised regions and the orbital 
angular velocity measured in "coordinate units" . These 
values for d and Q, together with the coordinate radius 
of the excision sphere r = 0.785798I37I, provide all pa- 
rameters necessary to reproduce the data in Table IV. 
The remaining quantities are dimensionless. mQo is the 
orbital angular velocity of the binary system as measured 
at infinity. E^/n is the dimensionless binding energy of 
the system with E\, defined as Eb = E ADM — m. J j \im 
is the dimensionless total ADM angular momentum of 
the binary system as measured at infinity, ijm is the di- 
mensionless proper separation between the two excision 
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TABLE IV: 


Sequence of corotating equal-mass 


black holes on 


circular orbits satisfying the Komar ansatz. 


The ISCO is at 


separation d 


~ 8.28. 














d 




rnQo 


Ei/fj, 


J/ fim 


i/m 




Sz/MC„ 




S K /M( TT 


40 


0.0052966 


0.01090 


-0.0233183 


4.9120 


22.87 


0.02148 


0.02150 


35 


0.0064212 


0.01327 


-0.0263265 


4.6620 


20.23 


0.02609 


0.02612 


30 


0.0080083 


0.01665 


-0.0301930 


4.4019 


17.58 


0.03264 


0.03267 


29 


0.0084050 


0.01750 


-0.0311008 


4.3488 


17.04 


0.03428 


0.03431 


28 


0.0088354 


0.01842 


-0.0320622 


4.2953 


16.51 


0.03606 


0.03610 


27 


0.0093039 


0.01943 


-0.0330816 


4.2414 


15.97 


0.03800 


0.03804 


26 


0.0098151 


0.02054 


-0.0341642 


4.1872 


15.43 


0.04012 


0.04017 


25 


0.010375 


0.02175 


-0.0353153 


4.1328 


14.89 


0.04245 


0.04250 


24 


0.010990 


0.02309 


-0.0365409 


4.0781 


14.35 


0.04501 


0.04507 


23 


0.011667 


0.02457 


-0.0378476 


4.0233 


13.81 


0.04785 


0.04791 


22 


0.012418 


0.02622 


-0.0392423 


3.9684 


13.26 


0.05099 


0.05106 


21 


0.013252 


0.02806 


-0.0407325 


3.9135 


12.72 


0.05449 


0.05457 


20 


0.014183 


0.03012 


-0.0423262 


3.8587 


12.17 


0.05841 


0.05850 


19 


0.015227 


0.03245 


-0.0440313 


3.8043 


11.62 


0.06282 


0.06293 


18 


0.016406 


0.03511 


-0.0458555 


3.7503 


11.06 


0.06782 


0.06794 


17 


0.017744 


0.03814 


-0.0478056 


3.6971 


10.50 


0.07352 


0.07366 


16 


0.019273 


0.04164 


-0.0498857 


3.6450 


9.943 


0.08007 


0.08024 


15 


0.021031 


0.04571 


-0.0520956 


3.5944 


9.377 


0.08765 


0.08785 


14.5 


0.022012 


0.04800 


-0.0532470 


3.5699 


9.092 


0.09190 


0.09212 


14 


0.023070 


0.05049 


-0.0544267 


3.5460 


8.807 


0.09650 


0.09674 


13.5 


0.024215 


0.05320 


-0.0556310 


3.5229 


8.520 


0.1015 


0.1018 


13 


0.025456 


0.05617 


-0.0568549 


3.5006 


8.231 


0.1070 


0.1073 


12.5 


0.026804 


0.05942 


-0.0580912 


3.4793 


7.941 


0.1129 


0.1132 


12 


0.028273 


0.06300 


-0.0593298 


3.4592 


7.649 


0.1195 


0.1198 


11.5 


0.029877 


0.06696 


-0.0605565 


3.4405 


7.356 


0.1266 


0.1271 


11 


0.031634 


0.07135 


-0.0617523 


3.4234 


7.060 


0.1346 


0.1350 


10.5 


0.033564 


0.07625 


-0.0628909 


3.4082 


6.762 


0.1434 


0.1439 


10 


0.035689 


0.08174 


-0.0639364 


3.3952 


6.462 


0.1532 


0.1538 


9.5 


0.038037 


0.08792 


-0.0648401 


3.3848 


6.159 


0.1642 


0.1649 


9 


0.040638 


0.09493 


-0.0655348 


3.3775 


5.853 


0.1765 


0.1773 


8.9 


0.041191 


0.09644 


-0.0656410 


3.3765 


5.792 


0.1792 


0.1800 


8.8 


0.041757 


0.09800 


-0.0657345 


3.3756 


5.731 


0.1819 


0.1828 


8.7 


0.042334 


0.09960 


-0.0658142 


3.3749 


5.669 


0.1847 


0.1856 


8.6 


0.042923 


0.1012 


-0.0658791 


3.3744 


5.607 


0.1876 


0.1885 


8.5 


0.043526 


0.1029 


-0.0659282 


3.3740 


5.545 


0.1905 


0.1915 


8.4 


fl 044 1 41 


fl 1 047 


-0.0659603 


"? 'WHS 






0.1945 


8.35 


0.044453 


0.1055 




3.3737 


5.451 


0.1951 


fl 1 Qfil 


8.3 


0.044769 


0.1064 


-0.0659741 


3.3738 


5.420 


0.1967 


0.1977 


8.25 


0.045089 


0.1074 


-0.0659738 


3.3738 


5.389 


0.1982 


0.1993 


8.2 


0.045412 


0.1083 


-0.0659685 


3.3739 


5.358 


0.1998 


0.2009 


8.1 


0.046068 


0.1102 


-0.0659419 


3.3743 


5.295 


0.2031 


0.2042 


8 


0.046738 


0.1121 


-0.0658929 


3.3748 


5.232 


0.2065 


0.2076 



surfaces as measured on the initial-data slice. 3 Finally, 
S z /M? r and Sk/M? t are two measures of the dimen- 
sionlcss spin of one of the black holes. S z is defined using 
the flat space Killing vector and Sk is defined using the 
approximate Killing vector. 



3 Note that the value of I listed in Tables IV and V of Ref.[5] 
were in error. The corrected values for the corotating case, and 
correct values for the "true" non-spinning case, are given in the 
following tables. 
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APPENDIX C: NON-SPINNING SEQUENCE 

In this Appendix, we list the numerical results for non- 
spinning equal-mass black holes in a quasi-circular or- 
bit as defined by the Komar ansatz. We have assumed 
conformal flatness, maximal slicing, and used Eq. (33b) 
for the lapse boundary condition on both excision sur- 
faces. The data is scaled relative to the total and reduced 
masses (m and /i) defined with respect to the irreducible 
mass of the apparent horizons M- m . 

In Table V, d and f^o are, respectively, the separa- 
tion of the centers of the excised regions and the or- 
bital angular velocity measured in "coordinate units" . 
f r = Q r /r2o is the rotation fraction necessary to obtain 
a non-spinning black hole. These values for d, Q, and f r , 



together with the coordinate radius of the excision sphere 
r = 0.7857981371, provide all parameters necessary to re- 
produce the data in Table V. The remaining quantities 
are dimensionless. mf2o is the orbital angular velocity of 
the binary system as measured at infinity. Eb/fJ, is the 
dimensionless binding energy of the system with Eb de- 
fined as Eb = -Eadm- m. J / frni is the dimensionless total 
ADM angular momentum of the binary system as mea- 
sured at infinity. Ijtn is the dimensionless proper sepa- 
ration between the two excision surfaces as measured on 
the initial-data slice. S Z /M? T is the the dimensionless 
spin of one of the black holes where S z is defined using 
the flat space Killing vector. The spin defined using the 
approximate Killing vector is zero by definition. 
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TABLE V: Sequence of non-spinning equal-mass black holes on circular orbits satisfying the Komar ansatz. The ISCO is at 
separation d ~ 7.55. 



d 


Oo 


fr 




Eb/fi 


J / lira 


t/m 




40 


0.0052973 


0.9551 


0.01090 


-0.0235329 


4.8717 


22.87 


-0.00004206 


35 


0.0064225 


0.9485 


0.01327 


-0.0266396 


4.6136 


20.24 


-0.00005974 


30 


0.0080108 


0.9396 


0.01665 


-0.0306751 


4.3423 


17.58 


-0.00008963 


29 


0.0084079 


0.9374 


0.01750 


-0.0316304 


4.2864 


17.05 


-0.00009813 


28 


0.0088388 


0.9351 


0.01843 


-0.0326458 


4.2298 


16.51 


-0.0001077 


27 


0.0093079 


0.9326 


0.01943 


-0.0337268 


4.1727 


15.98 


-0.0001187 


26 


0.0098198 


0.9299 


0.02054 


-0.0348800 


4.1150 


15.44 


-0.0001312 


25 


0.010380 


0.9270 


0.02175 


-0.0361122 


4.0568 


14.90 


-0.0001458 


24 


0.010996 


0.9238 


0.02309 


-0.0374317 


3.9979 


14.36 


-0.0001627 


23 


0.011676 


0.9203 


0.02457 


-0.0388473 


3.9386 


13.82 


-0.0001825 


22 


0.012428 


0.9165 


0.02622 


-0.0403693 


3.8786 


13.27 


-0.0002058 


21 


0.013264 


0.9123 


0.02806 


-0.0420092 


3.8182 


12.73 


-0.0002336 


20 


0.014198 


0.9077 


0.03013 


-0.0437799 


3.7574 


12.18 


-0.0002669 


19 


0.015247 


0.9025 


0.03246 


-0.0456957 


3.6962 


11.63 


-0.0003072 


18 


0.016431 


0.8967 


0.03512 


-0.0477726 


3.6348 


11.07 


-0.0003567 


17 


0.017775 


0.8903 


0.03816 


-0.0500280 


3.5732 


10.52 


-0.0004180 


16 


0.019313 


0.8829 
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